Coherent photon-hadron interactions in pA collisions: Small- a: physics after HERA 
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In this letter we study the photoproduction of heavy quarks and vector mesons in the coherent 
proton-nucleus (pA) interactions for RHIC and LHC energies and analyze if these processes can be 
used to determine the QCD dynamics at high energies. The integrated cross section and rapidity 
distribution are estimated using the Color Glass Condensate (CGC) formalism. A comparison with 
the linear dynamics predictions is also presented. Our results indicate that the nonlinear dynamics 
can be proven in those reactions, which are well suited for studing saturation effects. 



PACS numbers: 12.38.Bx; 13.60.Hb 



I. INTRODUCTION 

Over the last decade much progress has been realized 
towards understanding the QCD dynamics at high ener- 
gies. The successful operation of the DESY ep collider 
HERA has opened a new era of experimental and theo- 
retical investigation into the deep structure of the pro- 
ton and, in general, of hadronic interactions. Some of 
the most important observations are the striking rise of 
the proton structure function F2{x,Q^) for small values 
of the Bjorken variable x (< 10~^), the large contribu- 
tion of diffractive processes in this kinematical range and 
the geometric scaling (For a recent review, see e.g. Ref. 
0). Theoretically, at small x, due to the large gluon 
density, we expect the transition of the regime described 
by the linear dynamics, where only the parton emissions 
are considered, for a new regime where the physical pro- 
cess of recombination of partons become important in the 
parton cascade and the evolution is given by a nonlinear 
evolution equation. This regime is characterized by the 
limitation on the maximum phase-space parton density 
that can be reached in the hadron wavefunction (par- 
ton saturation), with the transition being specified by 
a typical scale, which is energy dependent and is called 
saturation scale Qsat (For recent reviews see Ref. Q)- 
Although the HERA experimental results have a natu- 
ral interpretation in terms of the saturation physics, due 
to the kinematical limitations of the experiment none of 
these phenomena can be taken as a conclusive evidence 
for a new regime of the QCD dynamics. Currently, the 
HERA II run has obtained much more precise experimen- 
tal data, but they are limited at center of mass energies 
smaller than 300 GeV, where deviations between linear 
and saturation predictions are small. Consequently, our 
understanding of the correct QCD dynamics at high en- 
ergies is still an open question. As HERA will stop its 
operations, in principle, in 2007, and the QCD dynamics 
must be known as precisely as possible in order to maxi- 
mize the discovery potential for new physics at the next 
generation of colliders, the study of alternatives which 



could constrain the QCD dynamics is timely and neces- 
sary. 

In this letter we study the photoproduction of heavy 
quarks and vector mesons in the coherent proton-nucleus 
{pA) interactions at RHIC and LHC energies and analyze 
if these processes can be used to determine the QCD 
dynamics at high energies (For similar studies on pp and 
AA collisions, see Refs. 0i 01)- The main advantage 
of using colliding hadron and nuclear beams for studying 
photon induced interactions is the high equivalent photon 
energies and luminosities that can be achieved at existing 
and future accelerators (For a recent discussion see Ref. 
i). 



II. COHERENT pA INTERACTIONS 

Lets consider the proton-nucleus interaction at large 
impact parameter [b > Rp + Ra) and at ultrarelativistic 
energies. In this regime we expect the electromagnetic 
interaction to be dominant. When we consider the elec- 
tromagnetic field associated to the proton and ion, we 
have that due to the coherent action of all protons in the 
nucleus, the electromagnetic field surrounding the ion is 
very larger than the proton one. This result can be easily 
understood if we use the Weiszacker- Williams method to 
calculate the equivalent flux of photons from a charge Z 
nucleus a distance h away, which is given by (For recent 
reviews see Ref. 3) 
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where lu is the photon energy, 7^ is the Lorentz boost 
of a single beam and 77 = Lob/^L', Koij]) and Ki{ri) are 
the modified Bessel functions. From the above expres- 
sion we have that photon spectrum of a nucleus with 
charge Z is proportional to Z^. Due to asymmetry in 
the collision, with the ion being likely the photon emit- 
ter, we have that the photon direction is known, which 
will implicate an asymmetry in the rapidity distribution 
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(see below). The coherence condition hmits the photon 
virtuaHty to very low values (Q^ < 1/R\), which implies 
that for most purposes, these can be considered as real. 
Therefore, coherent pA collisions can be used in order to 
study photon-proton interactions. 

The requirement that photoproduction is not accom- 
panied by hadronic interaction (ultraperipheral collision) 
can be done by restricting the impact parameter b to be 
larger than the sum of the proton and the nuclear ra- 
dius. Therefore, the total photon flux interacting with 
the target nucleus is given by Eq. Q integrated over 
the transverse area of the target for all impact param- 
eters subject to the constraint that the proton and the 
nucleus do not interact hadronically. An analytic ap- 
proximation for pA collisions can be obtained using as 
integration limit b > Rp + Ra, producing ^ 
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where = w {Rp + Ra)/il and Z^(?7) = (77) - (77). 
The cross section for the photoproduction of a final state 
X in a coherent pA collisions will be given by, 
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where ujmin = M^/i-fLnip, W^p = 2uo^/Skn and vSnn 
is the c.m.s energy of the proton-nucleus system. Con- 
sidering pPb {Ar) collisions at LHC, the Lorentz factor 
is 7l = 4690 (5000), giving the maximum c.m.s. en- 
ergy W^p w 1500 (2130) GeV. Therefore, while studies of 
photoproduction at HERA are limited to photon-proton 
center of mass energies of about 200 GeV, photon-proton 
interactions at LHC can reach one order of magnitude 
higher on energy. Consequently, studies of interac- 
tions at LHC could provide valuable information on the 
QCD dynamics at high energies. In this work we con- 
sider that the produced state X represents a QQ pair or 
a vector meson V . Since photon emission is coherent over 
the entire nucleus and the photon is colorless we expect 
that the events to be characterized by one (heavy quark 
production, with Y being the remaining of the proton) 
or two (vector meson production, with Y = p) rapidity 
gaps. 



III. QCD DYNAMICS AT HIGH ENERGIES 

The photon-hadron interaction at high energy (small 
x) is usually described in the infinite momentum frame 
of the hadron in terms of the scattering of the photon 
off a sea quark, which is typically emitted by the small- a; 
gluons in the proton. However, in order to disentangle 
the small-x dynamics of the hadron wavefunction, it is 
more adequate to consider the photon-hadron scattering 
in the dipole frame, in which most of the energy is car- 
ried by the hadron, while the photon has just enough 
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TABLE L The integrated cross section (event rates/month) 
for the photoproduction of heavy quarks and vector mesons in 
pA collisions at LHC (see text). 



energy to dissociate into a quark-antiquark pair before 
the scattering. In this representation the probing pro- 
jectile fluctuates into a quark-antiquark pair (a dipole) 
with transverse separation r long after the interaction, 
which then scatters off the proton 0. The main moti- 
vation to use this color dipole approach, is that it gives 
a simple unified picture of inclusive and diffractive pro- 
cesses. In particular, in this approach the heavy quark 
photoproduction cross section reads as. 



ct(7P- 
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where '^'[ Az, r) is the light-cone wavefunction of the 
photon 13 • The quark and antiquark helicities are la- 
beled by h and h. The variable r defines the relative 
transverse separation of the pair (dipole) and z (1 — z) is 
the longitudinal momentum fractions of the quark (anti- 
quark). The basic blocks are the photon wavefunction, 
vJ'T and the dipole-target cross section, adip- For pho- 
toproduction we have that longitudinal piece does not 
contribute, since I^'lP oc Q^, and the total cross sec- 
tion is computed introducing the appropriated mass and 
charge of the charm or bottom quark. 

Similarly, in the dipole picture the imaginary part of 
the amplitude for vector meson production at zero mo- 
mentum transfer reads as (See e.g. Refs. 0:I1|) 

ImA{-ip -> Fp) = ^ / dz d^r -^l^^ cTd^pii, r) ,(5) 
h,h 

where ^'J^^(z, r) is the light-cone wavefunction of the 
vector meson. The total cross section for vector meson 
photoproduction is given by 



cr (7P — + Vp) 



[ImA{s, t = 0)]^ 
I&ttBv 



(1 + /?') 



(6) 



where (3 is the ratio of real to imaginary part of the am- 
plitude and By labels the slope parameter. The values 
considered for the slope parameter are taken from the 
parameterization used in Ref. 0. 

We have that the total cross sections for vector me- 
son and heavy quark production in dipole approach are 
strongly dependent on the dipole-hadron cross section 
fdip, which contains all information about the target and 
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FIG. 1: (Color online) Rapidity distribution for heavy quark 
photoproduction on pA reactions for LHC energy (see text). 



the strong interaction physics. In the Color Glass Con- 
densate (CGC) formalism (Tdip can be com- 
puted in the eikonal approximation, resulting 



,r) = 2 



I [1 



S(x,r,6i)] 
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where S is the S'-matrix element which encodes all the in- 
formation about the hadronic scattering, and thus about 
the non-linear and quantum effects in the hadron wave 
function. The function S can be obtained by solv- 
ing an appropriate evolution equation in the rapidity 
y = ln(l/a;). The main properties of S are: (a) for the in- 
teraction of a small dipole (r ^ l/Qsat), S{r) « 1, which 
characterizes that this system is weakly interacting; (b) 
for a large dipole (r 1/Qsat), the system is strongly 
absorbed which implies S{r) <C 1. This property is as- 
sociate to the large density of saturated gluons in the 
hadron wave function. In our analysis we will consider 
the phenomenological saturation model proposed in Ref. 
[l2i| which encodes the main properties of the saturation 
approaches, with the dipole cross section parameterized 
as follows 



A/^(V) , forf<2, 

1 - exp [-a In^ (bf)] , for f > 2 , 



where f = rQsat(a;) and the expression for f > 2 (satura- 
tion region) has the correct functional form, as obtained 
from the theory of the Color Glass Condensate (CGC) 
0. For the color transparency region near saturation 
border (f < 2), the behavior is driven by the effective 
anomalous dimension jcS (a;, r) = jsa.t + '"iy^"* ; where 
7sat = 0.63 is the LO BFKL anomalous dimension at 
saturation limit. Hereafter, we label this model by CGC. 



IV. RESULTS 

The distribution on rapidity y of the produced final 
state can be directly computed from Eq. Q , by using its 
relation with the photon energy w, i.e. y cx \ii{2uj/mx)- 
Explicitly, the rapidity distribution is written down as. 



da[A + p^X + Y] dN^{uj) 



dy 



duj 



M-(8) 



Consequently, given the photon flux, the rapidity distri- 
bution is thus a direct measure of the photoproduction 
cross section for a given energy. In Fig. ^ we present 
our results for the heavy quark photoproduction at LHC 
energies. For comparison, we also present the predictions 
from the linear dynamics, denoted coUinear hereafter, 
which is calculated assuming that the coUinear factoriza- 
tion is valid and that the gluon distribution can be de- 
scribed bythe GRV98 parametrization (For more details 
see Ref. Q). In Tab. Uone presents the correspondent 
integrated cross sections (event rates), using a luminosity 
of £ppb = 7.4 X 10^^ cm^s~^. We have verified for com- 
pleteness that for RHIC energy the difference between the 
predictions for the rapidity distribution is small, which 
is expect due to small value of the photon-proton center 
of mass energy. However, at LHC we can observe a large 
difference between the predictions and having high rates. 
For instance, for charm quark CGC gives a cross section 
a factor 3 lower than coUinear and a factor 2 for bot- 
tom. This deviation holds even in case of experimental 
cuts on rapidity. Therefore, photoproduction of heavy 
quarks should provide a feasible and clear measurement 
of the underlying QCD dynamics at high energies. Since 
RHIC is obtaining data for dAu interactions and LHC is 
to be commissioned, these processes could be analyzed 
in the next years. The advantages are a clear final state 
(rapidity gap and low momenta particles) and no com- 
peting effect of dense nuclear environment if compared 
with hadroproduction. 

The present results can be directly compared with 
those obtained in Ref. which use only coUinear ap- 
proach. The orders of magnitude are similar for RHIC 
and LHC, with main deviations coming from different 
choices for the gluon pdfs and distinct factorization scale. 
Our estimates considering CGC (coUinear) for RHIC are 
142 (110) /ib for charm and 0.15 (0.10) /ib for bottom. 
The rapidity distributions can be also contrasted with 
estimations in Refs. |^|^, where the heavy quark pho- 
toproduction in AA and pp collisions have been com- 
puted. To see what is difference in the order of magni- 
tude among them let us perform a few parametric esti- 
mates. Roughly, the photon ffux on nuclei is approxi- 
mately Z"^ the flux on proton, dN^^/du cx Z'^dNP/dcu. 
Moreover, for heavy quarks have been not verified large 
nuclear shadowing ;3I such that (Tja ~ ^ cr^p • The ratio 
between pA production and pp (or AA) can be estimated 



as R 



■pA/pp(AA) 



da 



Therefore, using Eq. ||SJ), 
one obtains RpA/pp oc- Z^ and the enhancement reaches 
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FIG. 2: (Color online) Rapidity distribution for vector meson 
photoproduction on pA reactions at LHC energy (see text). 



FIG. 3: (Color online) Comparison between CGC and color 
transparency estimates (see text). 



a factor lO"*^ for heavy nuclei in the comparison between 
photoproduction on pA and in energetic protons. On the 
other hand, one obtains RpA/AA ^ ^-nd then pA 

is suppressed in relation to AA by a factor A. However, 
the larger pA luminosity, which is two order of magnitude 
higher than for AA, counteracts this suppression for the 
event rates. 

Concerning coherent meson production, in Fig. [21 one 
presents predictions for the rapidity distribution on pA 
collisions at LHC energy considering both light (Ar) and 
heavy (Pb) ions. These theoretical predictions are origi- 
nal in the literature. The corresponding integrated cross 
sections (event rates) are shown in Tab. having high 
rates. For completeness, we have computed the estima- 
tion for typical light meson (p) and a heavy one (J/*). 
For RHIC, the cross section are 13 mb (0.8 fih) for p 
( J/vt) production. These results can be contrasted with 
those obtained in the Refs. 0, 0, 0|. We can do a 
similar estimative of the ratio as for heavy quarks con- 
sidering that for light mesons a^A ~ A"^^^ o-jp and for 
heavy mesons a^A ~ A^^^ a^p Therefore, an en- 

hancement of for pA/pp remains but now RpA/AA 
yl-2/3 (yl-4/3-)^ respectively. For Pb ions at LHC this is 
a factor 35 for p and a factor 10"^ for J/\l/. It should be 
noticed that the coherent production of mesons is cur- 
rently measured at RHIC for the A A case. Therefore, 
the present estimation could be tested in dAu collisions 
with a good experimental feasibility. For comparison, 
we also present in Fig. [S] the prediction from the color 
transparency limit of the CGC model, which is obtained 
assuming that cr^^^ oc (rQsat)^ is given in all kinematic 
range by its value for f < 1. It is equivalent to disre- 
gard the saturation physics, which allows to estimate the 
importance of the CGC physics in the process. We have 
that the predictions for J/"^ production are somewhat 
similar, which is expected, since the heavy vector meson 



production is dominated by small pair separations, where 
the saturation physics does not contribute significantly. 
On the other hand, the photoproduction of p mesons is 
dominated by physics below saturation scale. The differ- 
ence is huge, reaching 3 orders of magnitude in rapidity 
distribution, which demonstrate the importance of the 
saturation physics on this process. 

Finally, lets comment on the uncertainties and model 
dependence for the results presented here. For the heavy 
quarks (HQ) production, the main theoretical uncertain- 
ties are the higher order corrections in perturbative ex- 
pansion and the quark masses. Within the perturbative 
QCD coUinear formalism, the higher order corrections 
for HQ photoproduction are quite small (a few percents) 
in contrast to the hadroproduction case and they can be 
taken into account by a suitable choice for the quark mass 
|l3lll5j |. The color dipole formalism, which we use in our 
analysis, contains contributions of higher orders by defi- 
nition as it is equivalent to the semihard factorization (or 
/ct-factorization) approach in leading logarithm approxi- 
mation. The latter is well known to include higher order 
perturbative QCD diagrams in its formulation. Concern- 
ing the DGLAP results (coUinear formalism), the uncer- 
tainties for different choices of the gluon distribution are 
small as it is probed at a relatively large scale p,"^ = 4mQ 
(approximately 9 GeV^ for charm and 81 GeV^ for bot- 
tom). Those uncertainties are known to take place at 
low virtualities = 1-3 GeV^. Moreover, our results 
are an extrapolation of the photoproduction cross sec- 
tion for energies beyond energies of current accelerator 
regime using the physical parameters which describe cor- 
rectly both the low energy (fixed target) and the HERA 
data points. Therefore, the present estimates should be 
reliable. Concerning vector meson production, in general 
it is described by quite distinct approaches for light and 
heavy mesons. In the light meson case, the old vector 
meson dominance (VMD) model is often adopted [l^ . 
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whereas for the heavy mesons a leading logarithmic ap- 
proximation of the coUinear approach is considered |l7j . 
In this case, the hard QCD scale is given by the vector 
meson mass. We have used the unique theoretical formal- 
ism available describing simultaneously light and heavy 
vector meson production Once again, the parame- 
ters of the model are fixed in order to describe the ex- 
perimental low energy and DESY-HERA datasets. The 
transition between light and heavy mesons is dynamically 
introduced by parton saturation effects (via saturation 
scale) in the proton target. 

V. SUMMARY 

The QCD dynamics at high energies is of utmost im- 
portance for building a realistic description oi pp/pA/AA 
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